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We use the diagram technique for spin operators to calculate Green's functions and observables 
of the spin- 1 quantum Heisenberg antiferromagnet on a square lattice. The first corrections to the 
self-energy and interaction are taken into account in the chain diagrams. The approximation repro- 
duces main results of Takahashi's modified spin-wave theory [Phys. Rev. B 40, 2494 (1989)] and 
is applicable in a wider temperature range. The energy per spin calculated in this approximation 
is in good agreement with the Monte Carlo and small-cluster exact-diagonalization calculations in 
the range < T < 1.2J where J is the exchange constant. For the static uniform susceptibility 
the agreement is good for T < 0.6J and becomes somewhat worse for higher temperatures. Never- 
theless the approximation is able to reproduce the maximum in the temperature dependence of the 
susceptibility near T — 0.9J. 
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I. INTRODUCTION 

Properties of the spin-^ quantum Heisenberg antifer- 
romagnet on a square lattice have attracted much at- 
tention in connection with the investigation of cuprate- 
perovskite high-temperature superconductors. Much of 
this interest stems from Anderson's original suggestion!!! 
that quantum spin fluctuation in CuC>2 planes of these 
compounds may be responsible for superconductivity. At 
present it is generally accepted that the undoped Cu02 
planes are well described by the Heisenberg model with 
nearest-neighbor interaction (see Ref. ^ and references 
therein). nn 

Numerical pCalculationsB'Q and the analysis of exper- 
imental datafl presented strong evidence that the two- 
dimensional nearest-neighbor s = h Heisenberg antifer- 
romagnet has long-range order at zero temperature. For 
T / the Hohenberg-Mermin- Wagner theoremQ shows 
rigorously that long-range order is destroyed. In real 
quasi-two-dimensional cuprate perovskites long-range or- 
der persists for temperatures lower than the Neel tem- 
perature. In this case the ordering is destroyed by small 
concentrations of carriers x < 0.01.0 In both cases there 
is a range of parameters where the arising short-range or- 
^ ■ der is characterized by the magnetic correlation length £ 
which is much larger than the in-plane intecsite distance 
a. In this situation, the spin-wave theorieaTLj modified 
for short-range order may be a go»d starting point for the 
investigation of antiferromagnetsJ13 For undoped crystals 
observables calculated in these theories agree nicely (cp. 
Refs. 0-||) with the small-cluster exact-diagonalization 
and Monte Carlo data up to temperature T ss 0.6 J where 
J is the exchange constaat. This corresponds to the cor- 
relation length £ w lOaE'S However, in cuprates for the 
region of carrier concentrations x •> 0.1, which is of prime 
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interest, £ is of the order of few lattice spacings.0 In this 
case the application of the modified spin-wave theories 
(MSWT) becomes doubtful. 

In this paper we extend the analytic description of el- 
ementary excitations of the undoped antiferromagnet to 
the region of higher temperatures and shorter correlation 
lengths. For this purpose we use the diagram technique 
for spin operators, developed in Ref. 12 (see also Ref . [l3[) . 
We consider the first corrections to the simplest chain di- 
agrams. Due to peculiar shapes of the main correction 
diagrams this approximation will be referred to as the 
one-loop approximation (OLA). We note that the OLA 
is not rotationally invariant (neither is the MSWT). 

As will be seen below, the obtained formulas are in 
many respects similar to the formulas of the MSWT de- 
veloped in Refs. f^,[| The major difference between these 
two approximations is in the evaluation of the excitation 
frequency. In the MSWT two parameters defining the 
frequency are deduced from the constraint of zero site 
magnetization and a self-consistency condition, while in 
the OLA these parameters are determined from the val- 
ues of diagrams. The OLA reproduces the main results 
of the MSWT and is applicable in a wider temperature 
range. The energy per spin calculated in the OLA agrees 
nicely with the Monte Carlo and small-cluster exact- 
diagonalization calculations up to T ~ 1.2 J. In the tem- 
perature range T < 0.6J both approximations give val- 
ues of the static uniform susceptibility in good agreement 
with numerical calculations. In the OLA, the agreement 
becomes somewhat worse for larger temperatures. Never- 
theless this approximation is able to reproduce the maxi- 
mum in the temperature dependence of the susceptibility 
near T = 0.9 J in close correspondence with experiment 
and numerical calculations, while in the MSWT the sus- 
ceptibility grows monotonously with temperature. 

It is worth noting also that the OLA reproduces re- 
sults of the traditional spin-wave approximation on the 
zero-temperature renormalization of the excitation, fre- 
quency and of the renormalization group theorju on 
the temperature variation of the correlation length. We 
have found also good agreement between data on zero- 
temperature spin correlations obtained with the pro- 
jected Monte Carlo simulations, the MSWT of Ref. |, 
and the present OLA. 

The outline of the paper is as follows. In Sec. II the di- 
agram technique for spin operators is discussed. The se- 
ries of the simplest chain diagrams is summed in Sec. III. 
At T — this approximation is similar tp, the spin- 
wave approximation suggested by AndersonJU The one- 
loop corrections to the chain diagrams are considered in 
Sec. IV. The obtained excitation spectrum is compared 
here and in Sec. V with the spectra of the traditional 
spin- wave theory and the MSWT. The spin- wave ap- 
proximation based on the obtained formulas is discussed 
in Sec. V. Here results of the OLA for the energy, suscep- 
tibility and spin correlations are compared with exact- 
diagonalization, Monte Carlo, and MSWT data. Our 
summary is given in Sec. VI. 
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II. DIAGRAM TECHNIQUE 



The Hamiltonian of the Heisenberg antiferromagnet 
can be written in the form 



H = JlmSlS„ 



lm 

Jim 

lm 



s \ s m + o ( s l^ s m + s l s m) 



(1) 



where 1 and m label the sites of the two sublattices of 
the square lattice, the exchange constants J\ m are ex- 
pected to be nonzero for nearest neighbor sites only, 

Jim = ^EaV+a witn a = (0, ±fl) , (±a, 0) , Si is tllC 

spin-i operator the components of which satisfy the com- 
mutation relations 

[«]=*4<W, (2) 

and analogously for cyclic permutations of indices x, y 
and z. Here n = 1 or m and s„ = sj± is v n . For s = j 
we have s„ = — i + s+s~ and s„s„ +s n 4 = 1. 

On one sublattice we change notations: s m = — s m 
and s m = — s m = i — s+s^. The new operators 
satisfy commutation relations (^|). Substituting the new 
notations in Hamiltonian (|l]) and omitting the tildes we 
find 



h = 2j fe^r+E 

\ 1 m 

H x = -£ J, 



• s i +s i + > v s+s m J , (.'-!) 

1 / j. j. . 



lm 



s t S l S m s m + r, l S l' S m + S l S mJ 



where N is the number of lattice sites. 

Our aim is the calculation of Green's functions 



D u ,(t) 


= -<PsT(t)s+) 


D[ m (r) 


= -(^ s r( r ) s m 


D'Ur) 


= -(Va+(T)a+ 



where sj~(r) = exp(/fr)s{~ exp(— Hr), the time-ordering 
operator V arranges operators in order of increasing time 
from right to left, and angular brackets denote thermo- 
dynamic averaging with the Hamiltonian H. It follows 
from the above definition that [-DL^t)]* = ^(r). 

Using Ho and H\ from Eq. (pto as the unperturbed 
Hamiltonian and perturbation, respectively, Green's fun- 
ctions can be represented by an infinite series generated 
by the known series of the evolution operator. The terms 
of the former series contain Green's functions 

(^(il'SW-'SMo ( 4 ) 

where the averaging and time dependences of the oper- 
ators are determined with the unperturbed Hamiltonian 
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Hq as indicated by the subscript of the angular brack- 
ets, so that 

As follows from Eq. (||), Green's functions (|J) contain 
only operators s + and (i.e. oti = + or — ) and do 
not contain operators s z . These Green's functions are 
nonzero when the number of s~ operators is equal to 
the number of s + operators. To calculate these functions 
we use Wick's theorem proposed for spin operators in 



Rcf. 12 (see also Ref. O). The theorem can be written 



in the form 

(VsZI(t 1 )...s-(t)...s%(t p )} 
= K(r - n )(V [s-,s%] ti . . . sZ(r p ))o + ■■■ 
+ K(t-t p )(Vs%(ti)...[s-, S %]J (5) 

where 

-2Jr f -i _ ,-2J/-A-! f 1, r > 



~ 2J P, T<0 



and = T is the temperature. The subscripts of 
the commutators in Eq. (^|) are the time arguments of the 
operators arising after the commutations. Equation (|5|) is 
easily verified by transferring s~ (r) to the right averaging 
bracket with the introduction of the commutators, using 
cyclic permutation in the averaging and transferring the 
operator from the left bracket to its initial position. 

As noted above, terms (|]) in the series expansion of 
Green's functions contain only operators s~ and s + . 
However, after the application of Wick's theorem (||) op- 
erators s z appear in the averaging. Equation (^|) is also 
applicable, if some of the s ai are s z operators and the 
couplings of with s z have to be taken into account 
along with the coupling of the former operators with s + . 
Equation ([5J) is used until only operators s z are left in 
the averaging brackets. Since the Hamiltonian Hq is the 
sum of terms related to individual sites, these averages 
are split into averages related to individual sites. Such 
averages are easily calculated using the recurrence rela- 
tion 



K»' + '^(«»')„K'.-^("y";„ 



~tanh(J/3) = -|. 



(6) 



For the following discussion we notice that 

( S n S n')o = ( S n)o — 2/3 <9J ^ S n)o ^nn'l (7a) 
( s n s n' s n")o = ( S n)o 

1 d 

~ ( S n>0 T^qJ (4>0 (<W + <W' + <W) 
1 d 2 



A(3 2 dJ 2 



(Sn) <5nn"W', (7b) 
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where we took into account that (s^) and its derivatives 
do not depend on n. Analogous formulas can be written 
for an average of a larger number of s z operators. 

The diagram technique can be— used to visualize the 
terms of the perturbation theory.113113 K(t — t') will be 
displayed graphically by a solid line with an arrow di- 
rected from r to t' [from the respective s~ to s + or s z 
operators, see Eq. (g)]. An s z operator is connected with 
the terminal point of this line unless a second line with an 
arrow enters into this point. This latter case corresponds 
to the situation when an s z operator which appears after 
the commutation of s~ and s + operators participates in 
another coupling. The s~ operator produced in this way 
participates in further couplings and therefore a third 
line has to issue from the terminal point of the former 



Excluding two free points labelled by the time argu- 
ments r and all other initial and terminal points of the 
directed lines are connected by wavy lines correspond- 
ing to the exchange integrals Ji m in the longitudinal [the 
first term in H%, Eq. (^)] and transverse interactions (the 
second and third terms in Hi). The transverse interac- 
tion line connects two directed lines, either terminating 
or initiating (Fig. [lb and c). The longitudinal interaction 
line connects four directed lines, one terminating and one 
initiating at each end of the interaction line (Fig. 0d). 

The delta functions in the last terms of Eqs. (pla) 
and (Qd) introduce limitations on the site indices of the 
s z operators. In the diagrams we shall denote these lim- 
itations by dashed lines connecting ends of directed lines 
(points corresponding to s z operators). A second-order 
diagram with such a dashed line arising due to the sec- 
ond term on the right-hand side of Eq. (|7|a) is shown in 
Fig. |l|e. As usual the partition function in the denomina- 
tors of Green's functions cancels all unlinked diagrams. 
It must be noted that a dashed line can connect an oth- 
erwise unlinked diagram, as shown in Fig. ^e. 



A satisfactory description of low-temperature excita- 
tions of the Heisenberg anti£crromagnet is achieved in 
the Anderson approximatiorJU which neglects the longi- 
tudinal interaction and substitutes the spin operators s + 
and s~ by Boson annihilation and creation operators. In 
the diagram representation this approximation is reduced 
to neglecting diagrams with longitudinal interaction lines 
and with dashed lines or loops formed by directed and 
transversal wavy lines (these two latter elements arise 
due to the non-Boson statistics of spin operators; it has 
to be noted that the loops contain the "triple points" 
similar to that shown in Fig. [l]a). The remaining chain 
diagrams for D\\i (r) are shown in Fig. ||. 

After the Fourier transformation 



two lines (see Fig. |T]a). 



III. CHAIN DIAGRAMS 
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where uj n — nTn is Matsubara's frequency with an even 
integer n, the series of chain diagrams can be written as 



D' c (k,iu; n ) 



where 

K(iuj n ) = - (iuj n - 2jy x , J k = 4J7 k 

are the Fourier transforms of K(t) and Ji m , respectively, 
7k = 5^ a exp(ika)/4 and the excitation frequency 

CJ k = 2J^1 - C 2 7 2. ( 9 ) 

At T — when c = 1 [see Eq. (p|)1 Eq. (^) reproduces 
the spin-wave spectrum of Ref. n4f For nonzero temper- 
ature c < 1 and a gap opens near k = 0. This gap is 
exponentially small for low temperatures J (3 1. 

IV. FIRST CORRECTIONS 

The first corrections to the chain diagrams contain 
one longitudinal interaction line, a loop with the "triple 
point" or a dashed line in the self-energy or in the renor- 
malized interaction. Similar corrections were considered 
for the case of a three-dimensional antiferromagnet and 
zero temperature in Ref. [l^ with the use of another dia- 
gram technique and for a ferromagnet in Refs. |T^ , p^| , ^6[ 
The respective diagrams are shown in Fig. [|. Diagrams 
(a) and (b) are self-energy corrections which are inserted 
in directed lines of the chain diagrams. Following Ref. [l5| 
we carried out partial summations in these diagrams: the 
line with two arrows indicates the series of chain diagrams 
(|J). Diagrams (d) to (f) are related to terminal points 
of directed lines and we include them together with dia- 
gram (c) to the renormalization of the interaction lines in 
the chain diagrams. However, in the chain diagrams the 
number of directed lines exceeds by one the number of 
wavy lines. Therefore the renormalized Green's functions 



= -cK{ico n ){ 1 + 



cJk 
2 



K(-iLu n )K(iuj n ) 



( 

V 2 

c (iujn + 2 J) 

(lUJ n ) 

- ^c 2 K (iuj n )K (~iuj n ) 



2 



K \—iu) n )K (iu>„ 



(8) 



2 

I c 2 J k 



I K(-iuj n )K(iuj n ) 



2 (ioj n ) 2 - Ull ' 



G 



have the additional multiplier a = 1 + erd + er e + £>/ where 
<7j is the contribution of the ith diagram in Fig. 0. Notice 
that the multiplier a was not included in Refs. ]l5] , |l6| . 
With these corrections Green's functions read 



ac (iuj n + 2 Je) 

£>'(k,iw n ) = 



D{\iw n) - — — -3 — — 
2 Jac 2 07k 



(10) 



(iw n ) - f2j 



where 



ft k = 2 J y e 2 - (c0 7k ) 2 , (11) 

£ = l + c 2 (l-/i)+c/, (12) 
^= l + c 2 J+^(l-di) 

-cosh" 2 (J/3)/' + J/3 cosh" 2 (J/3)/, (13) 
a = 1 + - c/i) - cosh" 2 (J/3)/' 

+J/3 cosh -2 (J/3)/, (14) 



/=-4£(-T-cx>th(^V (15a) 
A^4£J-c„ t h(^). (1 5b) 

Here and below summations over wave vectors are carried 
over the magnetic Brillouin zone. The second and third 
terms on the right-hand side of Eq. (|l^) are contribu- 
tions of diagrams (a) and (b) in Fig. ||, respectively. The 
second to fifth terms on the right-hand side of Eq. Jj"3] ) 
are contributions of diagrams (c) to (f), respectively. As 
mentioned, diagrams (d), (e), and (f) give the last three 
terms in a, Eq. ([hf). In the calculation of the correction 
introduced by diagram (c) we took into account that for 
any /(k) invariant with respect to operations of the crys- 
tal point group 

^k-k'Z(k') = JkJo -1 J k'/(k')- 

k' k' 

To estimate the contributions of the diagrams in Fig. |3| 
we introduce the function!] 

w ( x ) = Jjl>2 6 ( x -^)> 0<x<l, (16a) 

k 

and consider an infinite square lattice. In this case 
w(x) = JC (yi - x 2 ) = i + 0(1 - x), (16b) 
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where JC(k) is a complete elliptic integral of the first kind 
with modulus k. For low temperatures T <C J we find 



I = c- 2 (h-I 2 ) 



h 



w{x) 



: coth 



o yl — c?x 2 



w(l) + 2(1 - m ) 



1 — c 2 x 2 dx 



MX) 
pj 



ln(ij3J) + O [(/3jy 



(x)s/l - c 2 x 2 coth(/3J\/l - c 2 x 2 ^j dx (17) 
^(1)C(3) 



2(1 -mi) + 



2(/3J)3 
1 w(a;)a; 2 



■0[(/3J)" 



w(l)x 



In 



rdx 



1-Xq 



l + x 



where 



1 f uj( x\ 

= l-- / . V ; dx w 0.303398, 

2 Jo vl-i 2 

1 /•! 

1 = l-- / w(x)y / l-x 2 dx re 0.578974, 

2 Jo 



(18) 



x 



2Jc 



£(x) is the Riemann zeta function, and O [(/3J) ] in J' 
comprises terms slowly varying with oj n . 

From Eqs. @, (0), @, and (Jl§) one can see that 
due to the multipliers cosh -2 (J/3) for low temperatures, 
where the obtained formulas are expected to be valid, a e 
and <jf are much smaller than a c and ad and therefore the 
two former contributions can be omitted. The retained 
diagrams contain a closed loop formed by directed and 
wavy lines [diagrams (b), (c), and (d)] or a bubble [dia- 
gram (a)] . Due to these shapes of the main diagrams we 
use the term "one-loop approximation" . From the above 
equations we find that in two dimensions contributions 
of these separate diagrams are not small in comparison 
with the unperturbed self-energy and interaction. For 
example, the zero-temperature value of diagram (b) is 
2.38 J to be compared with the unperturbed self-energy 
2J. However, the sums of the values of the diagrams 
(a) and (b), as well as (c) and (d) are small in compari- 
son with the respective unperturbed values. Substituting 
( |l7| ) in Eqs. ([ll)-(|l4|) we find for an infinite lattice 



e = 2mi — 



2Je^l 



o [(pjy 



2{0Jf 

r = c^e -1 = 1 - 2m7/ 1 {l + 2w(l) (19) 

-4m + 2rrn + O [(PJ)- 1 ] } exp(-2 J(3), 
a = 2m a - w(l) + O [(pjy 1 ] . 
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The zero-temperature renormalization of the frequencies 
of elementary excitations, given by Eq. (|l9|), coincides 
with the result obtained previously in the traditional 
spin-wave approximatiorO£3 and by other methods. LtEI 
The temperature variation of the frequency described by 
the term -0.5ju(l)C(3)(/3 J)~ 3 in Eq. © is close to the 
MSWT resultO -0.322w(l)C(3)(/?J)- 3 . For nonzero tem- 
perature r < 1 which produces a gap near k = in the 
spectrum of elementary excitations. For low tempera- 
tures the size of the gap is exponentially small and thg 
obtained exponent — J/3 is close to the MSWT valueEl 
-4mo?7iiw" 1 (l)J/3 S3 -1.1037 J/3. However, the preex- 
poncntial factors are different in these two approaches. 

V. THE SPIN- WAVE APPROXIMATION 

In contrast to the small correction to the excitation 
frequency, the numerators of Green's functions dil]| ) are 
essentially renormalized in comparison with the chain ap- 
proximation (||). In Eq. ( |l9| ) a S3 —0.03 and with this 
small multiplier Green's functions ( [To| ) nearly vanish. It 
is worthj Hoting the direct analogy of this result with the 
MSWTQ'B where the correlations (sj~Sp), (sfs~) and 
the related Green's functions are exactly zero. Thus, like 
this theory, the OLA is not rotationally invariant and the 
z components of the spin correlations are much larger 
than the other components. In the following discussion 
these latter components will be neglected. 

The analogy with the MSWT may be continued, if 
attention is drawn to the fact that the quantities 

d(k, iuj n ) = D(k, iu n )/a, d'Qt, iuj n ) = D'(k, iuj n )/a 

are Green's functions of spin waves. Indeed, the fac- 
tor a in Green's functions ( |l0| ) differs from 1 due to 
the terminal-point diagrams (d), (e), and (f) which ac- 
count for the non-Boson statistics of spin operators cor- 
responding to the terminal points. Thus, dropping a in 
Green's functions (|l^) corresponds to the replacement of 
the spin operators by the respective Boson operators of 
spin waves. We introduce these latter operators with the 
Dyson-Maleev transformation 

sf=6 1 t , s+=(l-b\b l ) b h s{ = ^-b\b h 

s~ = -b s + = -tf (l — tfb ) s z = 1-6+ b 

(here and below we return to the initial coordinate sys- 
tem on the sublattice labelled by the index m). In 
these notations d and d' acquire the familiar form of the 
magnon Green's functions 

dw(r) = - {Vb\{r)b v ) , d[ m (r) = - (vb\(r)bl) . 

It is essential to note that the self-energy and interaction- 
line corrections, which are also connected with the non- 
Boson statistics of spin operators and renormalize the 







excitation frequency, are taken into account in d and dl . 
Analogous corrections are allowed for in the MSWT and, 
as a consequence, d and d' derived from Eq. (10) up to the 
factor c coincide in- their form with the magnon Green's 
functions obtainedtl in that theory. However, in the dis- 
cussed spin- wave approximation the excitation frequency 
is determined by Eqs._0)-(||), and @. They dif- 
fer from the equationsQffdcfiiiing the frequency in the 
MSWT. 

Taking into account that 



, \ f g ZT dz 

b\h) = -d lv (r = +0) = - I ——^d u ,(z) — 

Q) = - 



b v b{ 



-div(r 



e ZT dz 

—a Tdu'{z) — , 

e z P — 1 Zm 



where C is a closed contour embracing the complex plane, 
we find 



( SlS v) = {s! s f t )^(btbr)(b4) = 



c2 ^ 



oik(l'-l)_ 



2N 



coth J/kJl-rj^ 



-i 2 



(20) 



(sis m ) = (sfs z m ) « - (bib m ) = 



c 2 
2iV 



E elk(m " } "7=^ cothfj^yi - ^) 



In the derivation of Eq. (|20|) we have taken into account 
the constraint of zero site magnetization in the param- 
agnetic state 



(21) 



and in calculating (b^bnb^ba'J we have considered only 

diagrams corresponding to two independent magnon 
Green's functions, neglecting diagrams with interaction 
wavy and dashed lines between these two functions. 

The index k added to the parameter r in Eq. (|o|) is 
noteworthy. To fulfil constraint ( pi] ) we separate out the 
terms with k = 0, the so-called condensation parts, in 
the sums in spin correlations ( po| ) and suppose that r in 
the separated terms differs from the analogous parameter 
in terms with k ^ 0. This latter parameter is calculated 
from Eqs. @, @, @, and @, whereas r k=0 = r' 
in the condensation parts is determined from constraint 
(pl|) which can be written in the form 



iVyT 



,-'2 



: COth 



(jf3eV 



1 _ r >: 



'2 - 



1 + c 



(22) 



In a large lattice and at T — the condensation parts 
are equal to 2m , Eq. (|l^), and the sums in Eqs. j20| ) 
for large L — |L| are equal to 2m + 2 1 / 2 (irL)^ 1 . Thus, 
under these conditions 
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(s s L ) « 



m 



(23) 



where (— 1) L = +1 or —1 depending on whether the sites 
and L belong to the same or different sublattices. The 
value of the sublattice magnetization niQ w 0.3034 is 
good agreement with the Monte Carlo calculations 
The 1/L decay of the spin correlations to the square of 
the order parameter at zero temperature was expected in 
Rcf. [l^. For small L the zero-temperature spin correla- 
tions (|oSl) calculated from Eqs. (£§), @, ©, 
and (J22j), are compared with the results of the projected 
Monte Carlo methocD in Table |. The values agree nicely. 
We notice also that our spin correlations calculated to the 
fourth decimal place coincide with data of the MSWT of 
Rcf. |. 

For low temperatures the asymptotic behavior of the 
sums in Eqs. (Eof) for large distances is 



1 



d 2 k 



ikL 



TT 2 JpeJJ fc 2 + (2£)~ 2 
with the correlation length 

r 

t/a-. 



2 

J/3e 



:rL' M) { 2£ 



1 

4^2 V 1 + 2w (!) ^ 4m o + 2m i 



mi 



1/2 



,J0 



(24) 



Analogous exponential temperature dependences of the 
correlation length were obtained earlier in a number of 
works (see Refs. pp and references therein). The expo- 
nent J/3 in Eq. (£4) is close to the values 0.94 J (3 and 
27rmomi J (3 ~ 1.1037 J/3 obtained in Refs. || and |?], re- 
spectively. 

In Figs. |]^| the energy per spin 



E = 2J (s s a ) 
and the static uniform susceptibility 



(25) 



c 2 ^ 
12 



-£coth 2 ( J/3^1 



2 2 

?k7k 



(26) 



calculated with Eqs. @-((T|), @, ©, and 
compared with results of Monte Carlo simulations) 2 ^ 
the MSWT of-Jlefs. [?],|] and the exact diagonalization of 
a 4x4 lattice.0 The parameters r, e, and r' obtained in 
these calculations are given in Table ||. Respective pa- 
rameters rj and A which define the excitation frequency 
in the MSWTu w k = - 7? 2 7 2 [cf. Eq. ©] are also 
given in this table for comparison. The size dependence 
of x an d E calculated with the above formulas is negli- 
gible for large enough lattices — the difference between 
values obtained for a 64x64 lattice and those shown in 
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Figs. H and || for a 20x20 lattice is less than the size of 
the symbols in these figures. 

As seen from the figures, for temperatures T < 0.6J 
the results ab.tained with OLA are close to the results of 
the MSWTLra and are in good agreement with the Monte 
Carlo and exact-diagonalization data. From Figs. [| and 
^ we see also that energy ( p5| ) and the related short- 
range correlation are well described by the OLA even 
for T > J where the MSWT does not work. For the 
susceptibility in Fig. |^ the agreement is not so good for 
these temperatures, however the approximation is able 
to reproduce the maximum in x near T = 0.9 J and the 
decrease of the susceptibility for higher temperatures in 
close correspondence with experiment and Monte Carlo 
simulations. 

As seen from Table [ij, the parameters r' and e of the 
OLA are close to the respective parameters rj and A/ (2 J) 
of the MSWT of Ref. [?] for low temperatures. Essen- 
tial deviations begin at T « 0.5 J. We notice that the 
value of the parameter r is large at this temperature and 
grows with growing T . These values of r correspond to 
a large gap in the excitation spectrum near k = 0. To 
understand this result we notice that in the considered 
magnetic Brillouin zone the vicinity of the T point corre- 
sponds to two regions in the usual Brillouin zone which 
is twice as large as the magnetic one. These regions 
are located near the (0,0) and (n/a,n/a) points. Cor- 
respondingly, a branch of elementary excitations in the 
usual Brillouin zone transforms into two branches in the 
magnetic zone. In Ref. [| we have considered how the 
magnon branch, which is twofold degenerate at T = in 
the magnetic Brillouin zone, splits into two branches with 
growing temperature. These two branches form a unified 
branch in the usual Brillouin zone and are related to the 
central part and the periphery of this zone. In Ref. ^| we 
have indicated that the branch with the gap near k = 
in the magnetic Brillouin zone corresponds to the periph- 
ery of the usual zone and the states in this branch have 
larger spectral intensities in the magnon Green's function 
in comparison with states in the second, gapless branch. 
In the OLA and the MSWT of Refs. @§ the branch docs 
not split and for high temperatures its behavior is ex- 
pected to be determined by states with larger spectral 
intensities. Thus, the gap mentioned above is related to 
the vicinity of the (jr/a, it /a) point in the usual Brillouin 
zone. 



VI. SUMMARY 

In this.papjer we used the diagram technique for spin 
operatorsE3c3 to calculate Green's functions and observ- 
ables of the undoped spin-i quantum Heisenberg antifer- 
romagnet on a square lattice. We considered the first cor- 
rections - the one-loop diagrams - to the simplest chain 
diagrams. The obtained equations resemble the formulas 
of the modified spin- wave theory of Refs. ^,||. The ma- 
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jor difference between these two approximations is in the 
calculation of the excitation frequency — in the mod- 
ified spin-wave theory two parameters defining the fre- 
quency are deduced from the constraint of zero site mag- 
netization and a self-consistency condition, while in the 
one-loop approximation these parameters are determined 
from the values of the diagrams. The one-loop approxi- 
mation reproduces .the results of the traditional spin- wave 
approximatiorOlil on the zero-temperature renormaliza- 
tion of the escitation frequency, of the renormalization- 



spin correlations. Due to the mentioned difference in the 
excitation frequency and some other differences in formu- 
las, the one-loop approximation is applicable in a wider 
temperature range than the modified spin-wave theory. 
The energy per spin calculated in this approximation 
is in good agreement with the Monte Carlo and small- 
cluster cxact-diagonalization data up to the temperatures 
T ~ 1.2 J, whereas the modified spin-wave approximation 
gives quantitatively correct values of this energy in the 
range T < 0.6 J. In this temperature range both approxi- 
mations give values of the static uniform susceptibility in 
good agreement with numerical calculations. In the one- 
loop approximation, the agreement becomes somewhat 
worse for larger temperatures. Nevertheless this approx- 
imation is able to describe the maximum in the temper- 
ature dependence of this observable near T = 0.9 J in 
close correspondence with experiment and Monte_Carlo 
simulations, while in the modified spin- wave theoryu'B the 
susceptibility grows monotonously with temperature. 
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FIG. 1. Elements of diagrams. See text for an explanation 
of the different lines. 



FIG. 2. Chain diagrams for D u i(t) 



FIG. 3. The lowest-order corrections to self-energies and 
interaction lines. 



FIG. 4. rJjhe energy per spin obtained in the Monte Carlo 
simulation! 2 !] (•), in the modified spin- wave theory of Refs. |j|8| 
(o) and in the one-loop approximation of the present work 
(x). In the modified spin- wave and one-loop calculations a 
20x20 lattice was used. 



FIG. 5. The static uniform susceptibilita obtained in the 
Monte Carlo simulation for a 12 x 12 latticeEj (•), in the mod- 
ified spin-wave theory of Refs. (o) and in the one-loop 
approximation of the present work (x). In the modified 
spin- wave and one- loop calculations a 20x20 lattice was used. 



FIG. 6. The energ* per spin obtained by exact diagonaliza- 
tion of a 4x4 latticed (•) and in the one-loop approximation 
(x). 
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TABLE I. The zero-temperature spin correlations (s sl) 
obtained with the one- loop approximation (OLA) for a 20x20 
lattice m comparison with the projected Monte Carlo data 
(PMC)Jfl 



L 


PMC 


OLA 


(a,0) 


-0.3348 


-0.3354 


(a, a) 


0.2028 


0.2016 


(2a, 0) 


0.1772 


0.1751 


(2a, a) 


-0.1671 


-0.1648 


(2a, 2a) 


0.1475 


0.1454 


(3a, 0) 


-0.1491 


-0.1461 


(3a, a) 


0.1430 


0.1404 



TABLE II. Parameters r, e, r' in the one-loop approxi- 
mation (OLA) arid-the respective parameters of the modified 
spin- wave theoryD'B (MSWT). In both calculations a 20x20 
lattice was used. 





OLA 


MSWT | 


T/J 


1 - r 


e 




1 - 


r 




1 - 


- '1 


A/(2J) 


0.01 





8634E-04 


1 


1582 





4192E- 


04 





4193E- 


04 


1.1582 


0.1 





8633E-03 


1 


1583 





3188E- 


03 





3197E- 


03 


1.1583 


0.2 





1860E-02 


1 


1561 





6434E- 


03 





6480E- 


03 


1.1570 


0.3 





6692E-02 


1 


1490 





9904E- 


03 





1038E- 


02 


1.1516 


0.4 





2528E-01 


1 


1385 





1276E- 


02 





1566E- 


02 


1.1374 


0.5 





5953E-01 


1 


1257 





1508E- 


02 





2404E- 


02 


1.1097 


0.6 




0.1033 


1 


1101 





1735E- 


02 





4081E- 


02 


1.0629 


0.7 




0.1511 


1 


0928 





1976E- 


02 





8882E- 


02 


0.98897 


0.8 




0.1995 


1 


0752 





2230E- 


02 





2930E- 


01 


0.87371 


0.9 




0.2467 


1 


0589 





2498E- 


02 










1.0 




0.2918 


1 


0446 





2774E- 


02 










1.1 




0.3344 


1 


0327 





3055E- 


02 










1.2 




0.3743 


1 


0234 





3336E- 


02 
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Fig. 1. A. Sherman and M.Schreiber, One-loop approximation 




Fig. 2. A. Sherman and M.Schreiber, One-loop approximation 




Fig. 3. A. Sherman and M.Schreiber, One-loop approximation 




Fig. 4. A. Sherman and M.Schreiber, One-loop approximation 




Fig. 5. A. Sherman and M.Schreiber, One-loop approximation 




Fig. 6. A. Sherman and M.Schreiber, One-loop approximation 



